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HIGHLIGHTS 


•  A  new  method  is  proposed  to  measure  the  internal  temperature  of  (Li-ion)  batteries. 

•  The  temperature  measurement  is  based  on  electrochemical  impedance  spectroscopy. 

•  The  intercept  frequency  can  be  linked  to  the  internal  battery  temperature. 

•  No  hardware  temperature  sensors  are  required  and  heat  transfer  delay  is  absent. 
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A  new  method  is  proposed  to  measure  the  internal  temperature  of  (Li-ion)  batteries.  Based  on  elec¬ 
trochemical  impedance  spectroscopy  measurements,  an  intercept  frequency  (/0)  can  be  determined 
which  is  exclusively  related  to  the  internal  battery  temperature.  The  intercept  frequency  is  defined  as  the 
frequency  at  which  the  imaginary  part  of  the  impedance  is  zero  (Zjm  =  0),  i.e.  where  the  phase  shift 
between  the  battery  current  and  voltage  is  absent.  The  advantage  of  the  proposed  method  is  twofold:  (i) 
no  hardware  temperature  sensors  are  required  anymore  to  monitor  the  battery  temperature  and  (ii)  the 
method  does  not  suffer  from  heat  transfer  delays.  Mathematical  analysis  of  the  equivalent  electrical- 
circuit,  representing  the  battery  performance,  confirms  that  the  intercept  frequency  decreases  with 
rising  temperatures.  Impedance  measurements  on  rechargeable  Li-ion  cells  of  various  chemistries  were 
conducted  to  verify  the  proposed  method.  These  experiments  reveal  that  the  intercept  frequency  is 
clearly  dependent  on  the  temperature  and  does  not  depend  on  State-of-Charge  (SoC)  and  aging.  These 
impedance-based  sensorless  temperature  measurements  are  therefore  simple  and  convenient  for 
application  in  a  wide  range  of  stationary,  mobile  and  high-power  devices,  such  as  hybrid-  and  full  electric 
vehicles. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Li-ion  batteries  are  nowadays  the  technology  of  choice  to  supply 
portable  and  mobile  electronic  devices  with  electrical  energy.  This 
can  mainly  be  attributed  to  their  ability  to  combine  a  high  gravi¬ 
metric/volumetric  energy  and  power  density,  which  leads  to 
compact  and  low- weight  batteries  [1],  Since  the  90s,  Li-ion  batte¬ 
ries  have  come  to  dominate  small-scale  applications  such  as 
portable  phones  and  laptops.  However,  the  breakthrough  in  large- 
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scale  applications  such  as  (Hybrid)  Electrical  Vehicles  was  labo¬ 
rious,  which  is  mainly  caused  by  the  high  initial  cost  and  safety 
issues  [2],  To  guarantee  full  safety,  the  automotive  industry  has 
established  requirements  and  regulations  in  terms  of  security  and 
reliability.  Safety  for  batteries  (packs)  is  strongly  related  to  the 
temperature.  When  the  temperature  becomes  too  high,  the  battery 
can  catch  fire  or  even  explode.  For  the  success  of  the  electric  vehicle 
market  this  would  be  detrimental,  if  not  fatal  [3],  It  is  thus  of  high 
importance  to  have  a  reliable  on-board  thermal  management  sys¬ 
tem  with  proper  temperature  gauging. 

Surface  mounted  thermal  sensors  like  thermistors  or  thermo¬ 
couples  suffer  from  heat  transfer  delay  due  to  the  thermal  mass  of 
batteries  (packs).  Therefore,  if  the  system  is  not  in  thermal 
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equilibrium,  the  internal  battery  temperature  differs  from  that 
externally.  The  internal  temperature  can  forecast  the  risk  of  venting 
and  thermal  runaway,  which  might  lead  to  spontaneous  combus¬ 
tion  [4],  It  can  also  be  used  as  measurement  input  for  (thermal) 
Battery  Management  Systems.  Forgez  et  al.  [5]  measured  the  in¬ 
ternal  battery  temperature  using  a  thermocouple  to  validate  their 
model  and  to  determine  model  parameters.  A  thermocouple  was 
inserted  in  the  center  of  a  commercial  cell  by  drilling  a  hole  into  the 
battery.  This  is,  however,  not  a  suitable  approach  for  mass  pro¬ 
duction.  Lee  et  al.  [6]  developed  a  flexible  micro-temperature 
sensor  to  in-situ  monitor  the  temperature.  Obviously,  the  disad¬ 
vantage  of  this  method  is  that  a  piece  of  hardware  must  be  placed 
inside  the  battery. 

Methods  based  on  Electrochemical  Impedance  Spectroscopy 
(EIS)  do  not  suffer  from  these  disadvantages.  EIS  is  a  non¬ 
destructive  technique  designed  to  understand  the  operation  of 
electrochemical  systems  [7],  It  applies  a  harmonic  voltage/current 
signal  to  the  battery  and  analyzes  the  respective  current/voltage 
response  at  various  frequencies.  After  performing  EIS,  various 
diagnostic  plots  can  be  made.  In  particular,  Nyquist  plots  are  very 
instructive  and  may  reveal  a  sloping  line  at  low  frequencies,  a 
number  of  (depressed)  semicircles  at  mid  frequencies,  and  a  ver¬ 
tical  line  at  high  frequencies.  Although  a  lot  of  research  has  already 
been  conducted,  the  exact  identification  of  which  process  is 
responsible  for  a  specific  part  of  the  spectrum  is  not  trivial  [8-10],  It 
is  generally  accepted  that  the  high  frequency  semicircle  can  be 
attributed  to  the  surface  layer  on  the  active  electrode  material  [11  — 
18]  and  that  the  mid-frequency  semicircle  can  be  attributed  to  the 
electrical  double  layer  and  electrochemical  charge  transfer  reac¬ 
tion,  occurring  at  the  electrode/electrolyte  interface  [11,12,14-18], 
Sloping  lines  are  generally  attributed  to  diffusion  processes.  It  has 
been  reported  that  the  high  frequency  semicircle  of  commercial  Li- 
ion  cells  is  almost  independent  of  SoC  [9,12,14]  whereas  the  tem¬ 
perature  has  a  large  influence  [9,12,14], 

In  the  high  frequency  region  the  impedance  of  a  battery  is 
generally  dominated  by  its  inductive  behavior.  This  has  been 
attributed  to  the  porosity  of  the  electrodes  [19],  the  electrode  ge¬ 
ometry  of  spirally  wound  cells  [20]  or  the  conductive  path  formed 
by  the  terminals,  connectors  and  electrodes  [21—23],  As  the 
porosity  and  electrode  geometry  hardly  change  upon  cycling,  the 
inductance  is  expected  to  remain  more  or  less  constant  [10],  The 
inductance  of  commercial  Li-ion  batteries  was  found  to  be  inde¬ 
pendent  of  SoC  [14]  and  temperature  [9,14],  The  intercept  fre¬ 
quency  was  also  reported  to  be  independent  of  both  SoC  and  aging 
of  Li-ion  batteries  [24,25],  Although  it  is  known  from  literature  that 


the  intercept  frequency  changes  as  function  of  temperature  [9]  and 
is  independent  of  SoC  and  aging,  the  relation  between  these  has 
never  been  investigated  in  detail. 

In  this  paper  we  extend  on  our  earlier  work  [26],  in  which  we 
introduce  the  concept  of  gauging  the  temperature  by  measuring 
the  intercept  frequency.  The  intercept  frequency  is  defined  as  the 
frequency  at  which  the  Nyquist  plot  intersects  the  real  axis.  In  other 
words,  the  intercept  frequency  is  the  frequency  (Jo  [Hz])  or  angular 
frequency  (wo  =  2tt/o  [rad  s-1]),  at  which  the  phase  shift  between 
voltage  and  current  is  zero.  The  intercept  frequency  can  be 
measured  by  straightforward  four-point  impedance  measure¬ 
ments.  The  advantage  of  measuring  at  the  relatively  high  intercept 
frequencies  is  that  these  impedance  measurements  are  less  sensi¬ 
tive  to  variations  in  battery  load  currents.  This  is  due  to  the  fact  that 
the  intercept  frequency  is  beyond  the  frequency  range  where  the 
main  electrochemical  storage  reactions  of  batteries  take  place.  This 
makes  it  easier  to  perform  accurate  measurements  even  in  systems 
under  operation,  such  as  driving  electric  vehicles. 

2.  Theoretical  background 

A  schematic  representation  of  a  Li-ion  battery  is  shown  in 
Fig.  la.  Such  a  battery  consists  of  two  electrodes  which  are 
attached  to  a  current  collector  (Cu  and  Al  in  the  case  of  Li-ion)  and 
are  separated  by  a  separator  which  is  impregnated  by  a  Li-salt 
containing  electrolyte,  ft  is  generally  considered  that  a  Solid 
Electrolyte  Interface  (SEI)  is  formed  at  the  surface  of  the  negative 
and  a  surface  film  at  the  positive  electrode.  All  these  layers 
contribute  to  the  overall  battery  impedance.  A  comprehensive 
equivalent  electric  network  model  can  be  constructed,  as  shown  in 
Fig.  lb,  to  simulate  the  overall  impedance  where  L  represents  the 
inductances  and  Rs  the  ohmic  series  resistances  of  the  connec¬ 
tions,  electrodes,  electrolyte,  and  separator.  The  outer  parallel 
capacitors  and  resistors  represent  the  electrical  double  layer 
capacitance  (Cdi)  and  the  charge  transfer  resistance  (Ret)  of  the 
positive  and  negative  electrode,  respectively.  Two  more  RC- 
elements  represent  the  impedance  of  the  (SEI)  surface  layers. 
The  individual  capacitors  (Cd)  are  added  to  model  diffusion  in  both 
electrodes.  This  detailed  circuit  is  rather  complicated  for  an 
analytical  analysis  and  is,  therefore,  reduced  to  the  circuit  shown 
in  Fig.  lc.  In  this  reduced  circuit  a  single  RC-chain  is  combined 
with  an  inductance,  a  series  resistance  and  a  diffusion  capacitance. 
This  single  RC-chain,  comprising  R^n  and  Qdn,  describes  the 
combined  charge  transfer  kinetics,  double  layer  capacitances  and 
surface  layer  behavior  of  both  electrodes. 


L.H.J.  Raijmakers  et  al.  /  Journal  of  Power  Sources  247  (2014)  539-544 


541 


(dashed  lines)  mathematically  derived  from  Eq.  (5).  Insert:  theoretical  dependence  of 
fo  as  a  function  of  temperature.  Cd  =  1  F,  Ckm  =  0.01  F,  I  =  7  nH,  0  <  J?Mn  <  0.02  SI, 
A  =  400,000  Q-1,  Ej  =  10  kj  mol-1,  R  =  8.314  J  mor1  K  253  <  T  <  323  K. 


The  impedance  of  the  equivalent  circuit  shown  in  Fig.  lc  can 
easily  be  derived  from  the  impedances  of  the  separate  elements, 
according  to 

ZR  =  R,  Zc  =  -L  ZL  =  jaL.  (1) 

The  total  battery  impedance  (Zb)  is  then  given  by 


The  imaginary  part  of  Eq.  (2)  is  calculated  as 


(2) 


z.  uL  (Ckin  +  Cd)Ckin^in^  +  f 
Qf<^inKL'"2  •  l) 


(3) 


The  intercept  frequency  wo  can  be  calculated  by  equating  Eq.  (3) 
to  zero.  This  results  in  the  following  bi-quadratic  equation 


Q^o(Ck2iXn^O  +  l)  -  (Ckin  +  Q)CkilXnW2  _  , 

*=  CkinRkinCdiw0  +  [CdL  ~  (Qdn  +  QJQin^kin]  w0  ~  1  =  0 

(4) 

If  the  parameters  of  the  lumped  circuit  take  finite  positive  values 
then  Eq.  (4),  in  general,  has  four  (complex)  roots.  However,  since 
the  intercept  frequency  cannot  be  complex  and  must  have  a  posi¬ 
tive  value,  the  only  relevant  solution  of  Eq.  (4)  has  the  following 
form 


«(([a 


-(Cm  +  QXW$n1 

2^Xn  ] 

CdE-  (Ckin  +  Cd)Ckinf?gi|A 
2LCdCkinRkm  J 


Eq.  (5)  shows  that  wo  is  a  function  of  Rkin-  The  intercept  fre¬ 
quency,  fo  =  wo/(27t),  is  plotted  as  a  function  of  Rkin  in  Fig.  2./o  re¬ 
veals  a  sigmoid-shaped  curve  enclosed  between  two  limiting  lines. 
The  lower  and  upper  limit,  which  are  shown  as  dashed  lines  in 


Fig.  2,  are  mathematically  derived  from  Eq.  (5)  when  Rkin  ap¬ 
proaches  zero  and  infinity,  respectively.  It  is  generally  accepted  that 
Rkin  reveals  an  Arrhenius  behavior  [16,24],  which  can  be  given  by 


where  A  is  the  pre-exponential  factor,  £a  the  activation  energy,  R  the 
gas  constant  and  T  is  the  absolute  temperature.  When  Rkin  is 
substituted  in  Eq.  (5),  an  expression  for  f0  is  obtained  as  a  function 
of  temperature.  The  result  is  plotted  in  the  insert  of  Fig.  2  and  re¬ 
veals  that  a  decreasing  intercept  frequency  is  theoretically  expected 
at  higher  temperatures.  This  makes  the  intercept  frequency  suitable 
for  temperature  indication  of  the  battery.  In  the  following  sections 
this  theoretical  statement  will  be  checked  experimentally  and 
ready-to-use  algorithms  will  be  developed. 

3.  Experimental 

Three  batteries  have  been  investigated  experimentally:  a  pristine 
high  power  2.3  Ah  LiFePCU  ANR26650  cylindrical  battery  (A123 
Systems),  pristine  and  cycled  7.5  Ah  cylindrical  ultra-high  power 
Li(NCA)02  batteries  (GAIA  Akkumulatorenwerke  GMBH).  All  batte¬ 
ries  have  graphite  anodes.  All  pristine  batteries  have  been  activated 
by  (dis)charging  these  four  times  according  to  the  specifications  of 
the  manufacturers.  Charging  and  discharging  is  performed  using 
Maccor  automated  cycling  equipment  (Model  4200).  The  imped¬ 
ances  were  measured  in  the  galvanostatic  mode  with  an  Autolab 
PGSTAT30  (Metrohm)  in  the  frequency  range  of  0.01  Hz-30  kHz, 
using  50  frequencies  logarithmically  distributed.  The  current 
amplitude  was  set  to  200  mA  RMS.  During  all  measurements  the 
batteries  were  located  in  climate  chambers  (Maccor,  model  MTC- 
010).  Measurements  were  conducted  at  -20,  -10,  +10,  +30 
and  +50  °C.  Before  conducting  measurements  at  a  predefined 
temperature,  3  h  of  rest  have  been  applied  to  reach  thermal  equi¬ 
librium.  Five  SoC  values  have  been  investigated,  ranging  from  100% 
to  20%  in  steps  of  20%.  The  intercept  frequency  has  been  determined 
by  interpolating  the  imaginary  part  of  the  impedance  spectrum  to 

Zim  =  0. 

4.  Results  and  discussion 

Fig.  3  shows  the  Nyquist  impedance  spectra  of  a  pristine  NCA 
cell  at  -10  °C  and  varying  SoC.  The  total  impedance  varies 
considerably  with  SoC,  especially  the  large  semicircle  in  the  low 
frequency  region.  On  the  other  hand,  the  small  semicircle  and  the 
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Frequency  [Hz] 

Fig.  4.  Impedance  results  of  Fig.  3  represented  in  a  Bode  plot. 

inductive  tail  at  higher  frequencies  show  only  a  small  dependence 
on  SoC.  This  becomes  more  clear  in  Fig.  4  where  the  same  imped¬ 
ance  results  are  plotted  in  a  Bode  diagram.  The  intercept  frequency 
is  the  same  for  all  SoC. 

Fig.  5  shows  the  impedance  spectra  of  a  pristine  NCA  cell  at  60% 
SoC  at  various  temperatures.  As  expected,  the  large  semicircle  in¬ 
creases  significantly  at  lower  temperatures.  It  can  also  be  seen  that 
the  small  semicircle  exhibits  a  similar  behavior  as  the  large  semi¬ 
circle  and  that  the  whole  impedance  curve  is  shifted  to  the  left  as 
the  temperature  increases,  indicating  a  decreasing  ohmic  resis¬ 
tance.  Furthermore,  it  can  be  noticed  that  the  inductive  character¬ 
istics  remain  more  or  less  constant.  In  the  Bode  plot  of  Fig.  6  it  can 
be  seen  that  the  intercept  frequency  clearly  changes  as  a  function  of 
temperature.  It  is  evident  that  the  intercept  frequency  decreases 
with  increasing  temperature.  Plotting  the  intercept  frequency  as  a 
function  of  temperature  for  a  pristine  NCA  cell  (Fig.  7)  and  LiFePCU 
cell  (Fig.  8),  an  exponential  behavior  is  found  for  both  chemistries, 
which  is  essentially  independent  of  SoC.  These  measurement  re¬ 
sults  confirm  the  theory  described  in  Section  2.  The  intercept  fre¬ 
quencies  for  the  LiFeP04  cell  are,  however,  about  five  times  larger 
than  those  of  the  NCA  cell,  as  can  be  expected  for  different  sizes  and 
chemistries  (see  Eq.  (5)). 

From  these  results  it  can  be  concluded  that  a  similar  intercept 
frequency  dependence  on  the  temperature  is  found  to  be 


temperatures. 
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Fig.  9.  (a)  Voltage  discharge  curves  upon  cycling  a  NCA  cell;  (b)  Impedance  spectra  of  a  pristine  and  cycled  NCA  cell  at  60%  SoC  and  10  °C. 


Fig.  10.  Measured  Intercept  frequency  as  a  function  of  temperature  for  a  pristine  and 
cycled  NCA  cell  at  60%  SoC. 


independent  of  battery  chemistry,  dimension  and  storage  capacity. 
An  interesting  question  is  how  battery  aging  influences  the  inter¬ 
cept  frequency.  In  order  to  investigate  this,  a  7.5  Ah  NCA  cell  has 
been  cycled  for  about  600  times.  Some  of  the  voltage  discharge 
curves  are  shown  in  Fig.  9a.  About  25%  of  the  storage  capacity  is  lost 
and  the  discharge  overpotentials  are  significantly  increased  after 
600  cycles,  indicating  a  major  impedance  increase.  Fig.  9b  shows,  as 
an  example,  the  Nyquist  plots  of  the  pristine  and  cycled  NCA  cell  at 
10  °C  and  60%  SoC.  In  line  with  the  voltage  discharge  curves 
(Fig.  9a),  it  can  be  seen  that  the  impedance  of  the  cycled  cell  at  low 
frequencies  has  dramatically  increased.  In  this  particular  case  there 
is  no  significant  increase  in  series  resistance,  i.e.  increase  in  IR  drop. 
According  to  Eq.  (5)  an  increase  in  series  resistance,  however, 
should  not  influence  the  intercept  frequency  as  a  function  of  tem¬ 
perature,  because  it  is  not  present  in  this  equation  and  therefore 
has  no  influence  on  the  temperature  indication.  The  spectra  at  high 
frequencies  did  not  change  upon  cycling,  indicating  that  the 
intercept  frequency  is  independent  of  aging.  This  conclusion  is 
confirmed  in  Fig.  10  where  the  intercept  frequency  is  plotted  as  a 
function  of  temperature  for  both  the  pristine  and  the  cycled  cell. 
This  figure  proves  that  the  intercept  frequency  did  not  change  upon 
cycling. 


5.  Conclusions 


The  intercept  frequency,  i.e.  the  frequency  at  which  the  phase 
shift  between  voltage  and  current  is  zero,  is  extracted  from 
impedance  spectra  of  a  Li(NCA)02  and  a  LiFeP04  battery  and  related 
to  the  internal  battery  temperature.  A  theoretical  analysis  based  on 
a  simple  equivalent  electrical  network  model  revealed  that  the 
intercept  frequency  is  expected  to  decrease  with  increasing  tem¬ 
perature.  Experiments  confirmed  this  finding  and  additionally 
showed  that  the  intercept  frequency  does  not  depend  on  SoC. 
Measurements  on  a  cycled  NCA  battery  showed  that  the  intercept 
frequency  remains  stable  during  cycling.  Additional  measurements 
were  conducted  on  60  Ah  LiFePC>4  prismatic  batteries.  These  mea¬ 
surements  gave  results  similar  to  those  of  the  NCA  and  LiFePCU 
batteries.  The  proposed  temperature  indication  method  is  therefore 
readily  applicable  to  batteries  of  different  chemistries,  dimensions 
and  storage  capacities. 
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